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We study the dynamics of skyrmions under spin-transfer torque in Dzyaloshinskii-Moriya mate¬ 
rials with easy-axis anisotropy. In particular, we study the motion of a topological skyrmion with 
skyrmion number Q = 1 and a non-topological skyrmionium with Q = 0 using their linear mo¬ 
mentum, virial relations, and numerical simulations. The non-topological Q = 0 skyrmionium is 
accelerated in the direction of the current flow and it either reaches a steady state with constant 
velocity, or it is elongated to infinity. The steady-state velocity is given by a balance between current 
and dissipation and has an upper limit. In contrast, the topological <5 = 1 skyrmion converges to 
a steady-state with constant velocity at an angle to the current flow. When the spin current stops 
the Q = 1 skyrmion is spontaneously pinned whereas the <5 = 0 skyrmionium continues propaga¬ 
tion. Exact solutions for the propagating skyrmionium are identified as solutions of equations given 
numerically in a previous work. Further exact results for propagating skyrmions are given in the 
case of the pure exchange model. The traveling solutions provide arguments that a spin polarized 
current will cause rigid motion of a skyrmion or a skyrmionium. 


I. INTRODUCTION 

Soliton structures are found in ferromagnets and they 
can be considered as an encoding of magnetic informa¬ 
tion which is robust both under temperature and external 
probes. Stable topological solitons with the structure of 
a skyrmion had been predicted in the presence of the 
Dzyaloshinskii-Moriya (DM) interaction [lj [2] and they 
were observed in recent years as isolated structures m 
or forming lattices I5H23- In the presence of easy-axis 
anisotropy there are topological skyrmions with skyrmion 
number Q = 1 as well as non-topological Q = 0 solitons 
(27 t vortices) [H [9]. 

Skyrmions could be the stable and robust entities that 
are needed for the technology of recording and transfer¬ 
ring information, currently mainly obtained in magnetic 
media using domain walls [10] . The propagation of mag¬ 
netic information is done most conveniently by the injec¬ 
tion of electrical spin-polarized current. Single skyrmions 
and skyrmion lattices can be set in motion suggesting a 
promising technique for the manipulation of magnetic in¬ 
formation [EEHni- Propagation of skyrmions by spin cur¬ 
rent or the related spin-Hall effect may be a promising 
strategy for the implementation of racetrack memories 

m- 

The existence of two species of skyrmions (Q = 0 and 
Q/0) has allowed theoretical predictions for dramati¬ 
cally different dynamical behaviors pj|. We show here 
that spin torque accelerates a Q = 0 skyrmionium and we 
describe the process theoretically, while the study is com¬ 
plemented by numerical simulations. The skyrmionium 
may reach a steady state or it may absorb energy from 
the current and expand without limit. A skyrmionium 
propagating even when the external probe is switched off 
can be obtained. The situation is contrasted to the more 
well-studied case of a Q = 1 skyrmion under spin torque. 

The outline of the paper is the following. Sec. [n] 
gives a description of the Landau-Lifshitz model includ¬ 


ing damping and spin-transfer torques and provides the 
main theoretical tools. Sec. Ill is on the dynamics of a 
<5 = 0 skyrmionium and Sec. IV is on the dynamics of 
a Q = 1 skyrmion. Sec. [V] contains our concluding re¬ 
marks. An Appendix gives analytical results on the pure 
exchange model. 


II. MAGNETIZATION DYNAMICS UNDER 
SPIN-TRANSFER TORQUE 


We assume a thin ferromagnetic film with a 
Dzyaloshinksii-Moriya (DM) interaction. Let M(x,y,t) 
be the magnetization vector with M s the saturation 
magnetization and define the normalised magnetization 
m = M/M s , so that m? = 1. The conservative Landau- 
Lifshitz (LL) equation for the statics and dynamics of the 
magnetization is 

d t m = -m x / (1) 


and is valid in the absence of damping and external 
probes. We consider an effective field / which includes 
an exchange interaction with constant A , an easy-axis 
anisotropy perpendicular to the (aq, a; 2 )-plane of the 
film with constant K , and a DM interaction with con¬ 
stant I?P]. If the energy is W then the effective field 
/ = — SW/Sm is 


f = Am + k m 3 e 3 (2) 

- 2A [d 2 m 3 ei - d\ m 3 e 2 + {di m 2 — <9 2 toi) e 3 ] 


where we have used £d = 2A/|D| as the unit of length. 
The parameter 


K 


D 2 


K = 


= vw, Kn = — 


Kn 


4 A 


( 3 ) 


is the rationalized (dimensionless) anisotropy constant 
and A = D/\D\ = ±1 will be referred to as the chiral¬ 
ity. We choose chirality A = 1 in all of our numerical 
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calculations, while n is taken to be positive (easy-axis 
anisotropy). We have not included the demagnetizing 
field in Eq. 0 because it does not affect skyrmion con¬ 
figurations in a qualitatively significant way m ; it in¬ 
troduces a dependence of the skyrmion size on the film 
thickness [2T]. The time variable t in Eq. ([4]) is measured 
in units of r 0 = 2 AM s / (^D 2 ) where 7 is the gyromag- 
netic ratio. 

When a spin-polarized current is flowing in the plane 
of the film, say, in the X\ direction, and we also include 
damping effects, the magnetization obeys the Landau- 
Lifshitz-Gilbert equation with two additional terms for 
the spin-transfer torque 

(dt +udi)m = —m x f + m x ( ad t + f)ud\) m. (4) 



x, 


The dissipation constant is a while u is the effective spin 
velocity parallel to the spin current and f3 is the non- 
adiabatic spin-transfer torque parameter. 

In the absence of spin torque, that is, for u = 0 the 
ground state of the model is the spiral state for suffi¬ 
ciently small anisotropy. For n > k c = 7 r 2 /4 ss 2.4674 
the ground state is either of the two uniform ferromag¬ 
netic states m = (0, 0, ±1). In the latter case, skyrmions 
are excited states and they are classified by the skyrmion 
number defined as 

If 1 

Q=— qd 2 x, q = -e M „m • (0„m x <9 M m), (5) 

where q is called the topological density. The skyrmion 
number Q is integer-valued (Q = 0,±1,±2,...) for all 
magnetic configurations such that m = ( 0 , 0 , ± 1 ) at spa¬ 
tial infinity. For definiteness we will assume m = (0, 0,1) 
in all our calculations. 

Axially symmetric skyrmion configurations are con¬ 
veniently described in terms of the standard spherical 
parametrisation given by 

mi = sin© cos $, m 2 = sin 0 sin$, 771.3 = cos© ( 6 ) 

with the ansatz 

Q = e(p), $ = 0 + 71/2, (7) 

where (p, 0) are polar coordinates. Solving Eq. (nj) with 
boundary conditions 9(p = 0) = 7 r and 9(p —► 00 ) = 0, 
leads to a static skyrmion with Q = 1 shown in Fig. [T] If 
the boundary conditions are 9(p = 0) = 27r, 9(p —)• 00 ) = 
0 a Q = 0 configuration is found 0 which has been called 
a “skyrmionium” [19{[25] and is shown in Fig.[2j 

For the conservative LL Eq. 0 it can be proved that 
the moments of the topological density defined by 

4 = J x fJ/ qd 2 x, p = 1,2 (8) 

are conserved quantities EiHTj. The conservation laws 
hold in the case of an infinite film. For a magnetization 
field obeying an equation of the form d t m = — m x g the 
time derivative of the topological density is 

(9) 


FIG. 1: The axially symmetric (Q = 1) skyrmion represented 
through the projection (mi, m 2 ) of the magnetization vector 
on the plane. It is calculated as a static solution of Eq. 0 
for anisotropy k = 3. 



FIG. 2: The axially symmetric (Q = 0) skyrmionium repre¬ 
sented through the projection (mi, m 2 ) of the magnetization 
vector on the plane It is calculated as a static solution of 
Eq. § for anisotropy k = 3. 


as can be found by a straightforward calculation. For 
g = / it can be shown that = 0 and the ingredients 
of the proof for the specific / of Eq. (0) are given in 
Ref. H9j. 

Eq. Q including damping and spin torques can be 
written in the form 

dtm = — m x g ( 10 ) 

g = -—— y [f + am x f — (0 — a)ud\ m — a(/3 — a)um x dim] 
1 + or 

and the time derivatives for I^ can be found by using 
Eq. @. The moments are no longer conserved due to 


q = -e^d^(g ■ d„m) 
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the damping and spin torque terms and we find 

(1 + ot 2 )i\ = —(/? — a)u d \2 + a D 2 -T (1 + a(3)u ( 47 tQ) 
(1 + a 2 )i 2 = (/? — a)udn — aDi ( 11 ) 

where we have used the notation 


spin polarized current. In order to follow how the initial 
skyrmionium will be accelerated we will follow its lin¬ 
ear momentum which is defined via the conserved I^ of 
Eq. (|| as EHH1 

P^l = /i, v = 1 or 2. (17) 


d^ = J{d^m ■ d v m) d 2 x, 
Df, = J (m x /) • d^md 2 x. 


H,v = 1,2 


( 12 ) 


Eqs. © give an explicit result since they may be applied 
for any magnetic configuration. 

Let us consider a skyrmion which is initially static 
within the conservative LL Eq. 0. and we suddenly apply 
an electrical current according to Eq. 0 . This will start 
to move and the overall motion is given by Eqs. (11). As 
a next step we will assume that it will eventually reach 
a steady-state with velocity v = (vi,v 2 ). We may then 
write the traveling wave ansatz for the magnetization 


m{x 1 ,x 2 ,t) = m 0 (£i, 6 ;ui,v 2 ) (13) 

= Xi - V!t, £, 2 =x 2 - v 2 t 


so that dtm = —v v d v m with v = 1, 2. Inserting this in 
Eq. 0 we obtain 


ud\m — v v d v m = m x f + m x (/3 uditn — av u d u m). 

(14) 

We will assume that d \, d 2 in the above equation denote 
derivatives with respect to £ 1 ,^ 2 - 

We may now take the cross product of both sides in 
Eq. (14) with d^m for /i = 1 or 2 and then contract 
with m. In the result, the term containing the effective 
field / (i.e., the term due to the conservative part of the 
equation) is written as a total derivative ( 2 U EZ] : 


- f ■ d^m = d u a^ 


(15) 


where the explicit form of the tensor for the effective 
field 0 is given in Ref. jT5| . Upon integrating over all 
space the total derivative ( |15[ ) vanishes and we obtain the 
pair of virial relations [55] ns 


(— 47 tQ + ad 2 i)vi + ad 22 v 2 = fiud 2 \ 
aduvi + ( 47 tQ + ad\ 2 )v 2 = fhidn. 


u ( 47 rQ) 


(16) 


Some general conclusions can be drawn from Eqs. © 
and |l 6 | ). The dynamics and the steady-state velocity 
(t'lj u^Jcan be obtained in particular cases as will be dis¬ 
cussed in the next sections for the cases of a skyrmionium 
and a skyrmion. 


III. A TRAVELING Q = 0 SKYRMIONIUM 

Let us consider the Q = 0 skyrmionium of Fig. [2] which 
is a static solution of Eq. (|T]) and we suddenly apply a 


To be sure, the above quantities have the meaning of a 
linear momentum within the hamiltonian equations ( 1 ). 
We will though extend their use in the full model (4). 
The time derivatives of the components of the linear mo¬ 
mentum are given by Eqs. (|ll|) applied for Q = 0: 


(18) 


(1 + a 2 )P\ = (P — a)u dn -«5i 

(1 + a 2 )P 2 = ((3 - a)u di 2 - aD 2 . 


These relations give a simple result when we apply them 
for the initial skyrmionium which is a static solution of 
Eq. Q (thus D^ = 0 ) and is axially symmetric (thus 
di 2 = 0). We obtain 


Pi = 


— a 


u di 


P 2 = 0 . 


(19) 


The skyrmionium is accelerated acquiring a linear mo¬ 
mentum component along the Xi-direction only. The ac¬ 
celeration is zero when (3 = a and this point will be 
clarified in the following when we present solutions of 

Eq. Q. 

If the skyrmionium eventually reaches a traveling 
steady state then the virial relations (16) apply. For 


0 = 0 they reduce to a simple form and give the velocity 
of the steady state: 


dn di 2 \ f avi - /3u\ = (0 
d 2 i d 22 ) l av 2 J 10 



( 20 ) 


provided det(ci Al!/ ) ^ 0. Therefore, the skyrmionium in 
steady state moves in the direction of the current. The 
presence of dissipation is crucial as there is apparently no 
steady state for a = 0. The acceleration process along 
the axis of the current described by Eqs. (191 is compat¬ 
ible with the steady-state velocity in Eq. (20 1 . 


We can find exact traveling solutions for a skyrmion¬ 
ium under spin-transfer torque. In order to show how 
this can be achieved, we start by considering the special 
case /3 = a, for which Eq. Q becomes 


)m = — 


m x f + am x (d t + udi) m. ( 21 ) 


We look for traveling wave solutions of the form ([13j) and 
we choose Vi = u, v 2 = 0. We have d t m = —udim, 


which is used to reduce Eq. (21) to m x / = 0. Thus, 
if we choose n(xi,x 2 ) to be a static solution of Eq. 0 
then m(xi,x 2 ,t) = n(£i, X 2 ) gives a configuration which 
satisfies Eq. (21) and is a traveling solution with velocity 
(vi,v 2 ) = (m, 0). In conclusion, the static skyrmionium 
solution of the conservative LL Eq. § (shown in Fig. 0 
is a traveling solution of the full equation 0 for the case 
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FIG. 3: The velocity components (vi,v 2 ) = ( Xi,X 2 ) for 
a skyrmionium under spin torque with parameter values 
( |25| ). The expected velocity at the steady state is (y\,v 2 ) = 
(0167,0). 


0 = a. That also explains the vanishing acceleration in 
Eqs. (19). This mathematical result has the following 


physical content. If we apply spin polarized current to 
an initially static skyrmionium this is expected to be set 
into rigid motion without significant deformations. The 
same can also be argued for the traveling skyrmion in the 
next Section. 

Let us now generalize the above for 0 ^ a. The trav¬ 
eling wave ansatz ( |13|) with the choice (vi,v 2 ) = (u, 0) is 
used to reduce Eq. Q to 

(zt — v)d\m = —m x / + (0u — av) m x dim. ( 22 ) 

We now choose v = flu/a and the equation is further 
reduced to 


vo dim = m x f, 


vo = 


(23) 


This equation is identical to that for a steady-state travel¬ 
ing with a velocity u 0 within the conservative LL Eq. (JT]) . 
Such states were numerically calculated and studied in 
Ref. 1191 . A family of traveling skyrmioniums were found 
with velocities up to a critical velocity v c ss 0.102. In 
conclusion, if we denote the solutions of Eq. ( [23] ) by 
n(xi,X 2 ',Vo) then the form m(xi,x 2 ,t) = n(£i,x 2 ;zV) 
with = xi — vt is a traveling wave solution of Eq. Q 
with velocity v = 0u/a. The condition vq < v c for the 
skyrmionium configuration satisfying Eq. (23), becomes 
in the presence of spin-transfer torque 


v < u + v c 


(24) 


We have conducted a numerical simulation using the 
parameter set 



FIG. 4: Contour plots of m 3 for a skyrmionium under spin 
torque with parameter values (251. (Left) The initial condi¬ 
tion, at t = 0, is a static skyrmionium solution of Eq. 0 - 
(Right) The skyrmionium at t = 40 when it has been accel¬ 
erated. The inner part has moved down relative to the outer 
part. The contour levels plotted are m 3 = 0.9, 0.6,0.3,0.0 
(solid lines) and m 3 = —0.3, —0.6, —0.9 (dashed lines). 


We use as an initial condition the skyrmionium of Fig. [2] 
with k = 3 and apply the spin current. The expected 
velocity in the steady-state is given by Eq. (20) and is 
= 0.167, 


Vi = U.ibC, v 2 = 0. The prediction is confirmed by 
the numerical simulation. We define the position of the 
skyrmionium as (Xi,X 2 ) with 


/ x /J (l - to 3 ) d 2 x 


M f( 1 — m 3 ) d 2 x 


(26) 


Fig.[3]shows the velocity (zq, v 2 ) = (Xi,X 2 ) as a function 
of time. The skyrmionium is accelerated and the change 
of linear momentum at t = 0 verifies the prediction P\ = 


0.34 calculated by Eq. (19) when we use the numerically 
calculated value du = 57.1. The velocity at t = 0 is 
vi(t = 0 ) « u = 0.1 and it increases to v\ (t = 100 ) = 0.16 
at the end of this simulation. The component v 2 acquires 
some small value at the initial stages of the simulation 
and it later goes to zero. The results have been confirmed 
also by a simulation in a moving frame, running for times 
longer than those shown Fig. [3] where it is seen that 
(zq,z; 2 ) converge to the expected values at the steady- 
state. 

Fig. a shows two snapshots of the skyrmionium under 
spin torque. At the initial time t = 0 we have the axi¬ 
ally symmetric static solution of Eq. ([l]). The accelerated 
skyrmionium at t = 40 has velocity zq = 0.143 and has 
lost axial symmetry: its central part has moved lower. It 
is very similar to the propagating skyrmionium studied 
in Ref. (19]. We conclude that the application of spin cur¬ 
rent is a method to obtain a propagating skyrmionium in 
a steady state. We note that the skyrmionium continues 
to travel at its acquired velocity when the spin currents is 
switched off, irrespectively of whether a steady state was 
reached or not. This is in stark contrast to the dynamics 
of a skyrmion or to ordinary domain wall dynamics. 

Let us now consider a second set of parameter values 


(25) 


a = 0.04, 


(27) 


a = 0.06, 


u = 0 . 1 , 0 = 0.1. 


u = 0 . 1 , 0 = 0.1. 
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which gives a velocity for the skyrmionium v = f3u/a = 
0.25 violating condition (24). In this case the skyrmion¬ 
ium is accelarated until its velocity approaches the limit¬ 
ing value u+v c while the configuration becomes elongated 
along the X 2 axis and eventually reaches the boundaries of 
our numerical grid. Presumably, the process would con¬ 
tinue until the skyrmionium configuration! is destroyed 
or until it turns to a domain wall extending to infinity in 
the X 2 direction. 


IV. A TRAVELING Q = 1 SKYRMION 

Let us now consider topologically nontrivial solutions 
(Q 7 ^ 0) such as the Q = 1 skyrmion of Fig. [T] which is 
a static solution of Eq. ([I]), and suddenly apply a spin 
polarized current. In order to follow the skyrmion as it 
moves we will follow the coordinates of its guiding center 
defined as the normalized moments in Eq. ([8]): 

= 4 (28) 

They give a measure of the position of a skyrmion and 
are conserved quantities as explained in connection with 
Eqs. 0, (0. 

The instantaneous velocity (i?i,I?2) is given through 
Eqs. © For the initial axially symmetric skyrmion 
solution of Eq. 0) we have D^ = 0 and di2 = 0, d\\ = 
d ,22 = W ex , where W ex is the exchange energy. Eq. © 
gives 


/3 = a since it emerges again as a special case as seen in 
Eq. (29). In a steady state Eq. (30) gives for the skyrmion 


a velocity (iq, V 2 ) = (u , 0) colinear with the spin current. 
When we substitute this in Eq. (14) we obtain the static 


LL equation mx f = 0. Thus the argument employed in 
Sec. |III| for a skyrmionium also applies for a skyrmion: a 
static skyrmion solution of the conservative LL Eq. (jTj) , 
which we denote n(x 1,22), is a traveling solution of the 
full equation Q with m(xi,X 2 ,t) = n(£ 1,2:2) and £1 = 
x\ — ut. 

It is well-known that there are no traveling skyrmion 
solutions of E q. (|Tj ). i.e., there are no Q / 0 skyrmion so¬ 
lutions of Eq. ( |23| ), and this can be rigorously established 
[26]. Therefore, the arguments about traveling skyrmio¬ 
nium solutions following Eq. (23) cannot be applied in 
the case of skyrmions. 

For ^ / a a traveling skyrmion should have 1)2 / 0 
as shown by Eq. (30), that is, the skyrmion travels at 
an angle with respect to the direction of the flow of the 
spin current. In the case /3 7^ a we could not find exact 
traveling wave solutions of Eq. Q , for the effective field 
©■ Exact results are though indeed obtained for the 
pure exchange model in Appendix [A] 

For small deviations from the simple case (3 = a and 
(vi,V 2 ) = (u,0) we may assume that the skyrmion re¬ 
tains approximately axial symmetry, and thus we have 
a diagonal d = dS^ where d is a constant. Eqs. (30) 
have now a relatively simple solution: 


V\ = u 


1 + a/3d 2 
1 + (ad) 2 


= u + ad V 2 , V 2 = u 


(/3 — a) d 
1 + (ad) 2 


(32) 


6 l + a/3 , a h 

Ri = u ——= u + -=R 2 ■ 
1 + a 2 d 




where we denoted d = dn/(47r<2) = d22/(47rQ). Thus, 
the skyrmion will initially have a velocity component in 
the direction of the current flow while its velocity com¬ 
ponent perpendicular to it depends on the sign of /3 — a. 

If we now assume that the skyrmion will eventually 
reach a propagating steady state then its velocity will 
satisfy the virial relations (161. For Q / 0 we define 


d^ v = d^vI(dmQ) and write the virial relations as 


For the pure exchange model d = 1 (see Appendix 0 
while for other models such as in Eq. ([2]) we have d > 1 
calculated by substituting the static skyrmion solution of 
Eq. Q in Eq. (12). Eq. (32) gives the so-called mobility 
relation, i.e., a linear relation between the velocity and 
the current. 

In order to check the theoretical predictions we have 
conducted a numerical simulation using the parameter 
set 


a = 0.2, u = 0.1, 


= 0.5 


(33) 


(1 — ad\ 2 )vi — ad 22 'V 2 = (1 — (dd\ 2 )u 
aduvi + (1 + ad\ 2 )v 2 = /3udn. 


(30) 


These imply that, in general, both components of the 
velocity are nonzero, vi,V 2 7^ 0, therefore a propagating 
skyrmion will move at an angle to the flow of the spin 
current. Unlike in the case of a skyrmionium, dissipation 
is not necessary in order to obtain a steady-state in the 
case of the skyrmion and Eqs. (30) give for a = 0 the 
velocity 


Vi=u — fiudi2, V2 = fiud\i. (31) 

In order to obtain more detailed information on the 
propagating skyrmion configuration we study the case 


with large values for a, fj and /3 — a. We use as an ini¬ 
tial condition the skyrmion of Fig. [l] with k = 3 and 
apply the spin current. The numerically found veloc¬ 
ity for the skyrmion is shown in Fig. [5] We initially 
have R\ = 0.106, R2 = 0.044; the velocity presents 
oscillations and eventually converges to constant values 
(^1,^2) = (0.1140,0.0436). The expected initial velocity 
is found from Eqs. p9| ) and it is R\ = 0.106, R2 = 0.046, 
where we have used the numerically calculated value for 
the exchange energy W ex = 20.08 => d = 1.598. The 
final velocity is in excellent agreement with the veloc¬ 
ity at a steady state predicted by Eq. (32) which gives 
(ui, U2) = (0.1139, 0.0435), where we assume that the ini¬ 
tial skyrmion profile is not significantly changed. The 
skyrmion configuration is indeed not visibly distorted 
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On the other hand, the skyrmionium motion presents 
some similarities with the skyrmion motion in a stripe 
geometry m- 

Our methods can be also applied to simpler one- 
dimensional models (wires) provided static domain wall 
solutions of Eq. ([l]) exist. 


Appendix A: The pure exchange model 

If we set f = Am we have the so-called pure exchange 
model and we will present analytic results which elu¬ 
cidate the discussion in Sec. [TV In the pure exchange 
model the Bogomolnyi relations Unim] 


FIG. 5: Simulation results for the velocity components 
(ui, V 2 ) = Ri,Ri) as a function of time t for a skyrmion under 
spin torque with parameter values 


during the simulation compared to the initial axially sym¬ 
metric skyrmion. 

We have also conducted a numerical simulation using 
the parameter set (27) and the results are again in ex¬ 


cellent agreement with the predictions of Eqs. (29) and 

(321. 


We finally note that the skyrmion is pinned at its final 
position when the spin current is switched off. This is 
contrasted with the dynamics of a skyrmionium which 
travels at a constant velocity also in the absence of ex¬ 
ternal forces (and damping). 


V. CONCLUDING REMARKS 


d\m = m x 82 m, 82 m = — m x 8 ±m, (Al) 

which contain only first order derivatives, are sufficient 
for obtaining static solutions of Eq. ([I]), i.e., solutions for 
m x / = 0. A large class of Q yf 0 skyrmion solutions can 
be found by solving (Al). Of those, the axially symmetric 
<5=1 skyrmion configuration of the form 0 will be 
denoted m = n( xi,x 2 ) and reads 


n 1 = — 


2 CIX 2 


n 2 = 


2ax\ 


n 3 = 4t4> (A2) 


where a is a arbitrary positive constant giving the 
skyrmion radius and p 2 = xf + x%- Configuration (A2) is 


similar in its gross features to that shown in Fig. [T]"" 

We turn to Eq. (I4| ) for a traveling steady state with 
velocity (t’i, V 2 ) and we require 


u — v 1 = — av 2 
V 2 = fiu — av\ 


V! =^U 
, V2 = 1+R5 u ' 


(A3) 


We have studied the dynamics under spin torque 
of non-topological (Q = 0) and topological (Q yf 0) 
skyrmions in films of Dzyaloshinskii-Moriya materials 
with easy-axis anisotropy. Analytical results are obtained 
using the equations for the linear momentum (for Q = 0) 
or the guiding centre (for Q / 0), and virial relations. 
The study is complemented by a set of numerical sim¬ 
ulations. Furthermore, we obtain exact solutions of the 
Landau-Lifshitz equation including spin torques, given in 
Eq. Q, for some particular cases. The traveling solutions 
obtained provide arguments that a spin polarized current 
will cause rigid motion of a skyrmion or a skyrmionium. 
This result is applicable not only for the solitons studied 
in the present work but also for other models, e.g., for 
driven motion of magnetic bubbles. 

The Q = 0 skyrmionium is accelerated by the spin 
torque and it continues moving after switching off the 
current. This Newtonian dynamics was also observed 
in the case of the skyrmionium under an external field 
gradient m ■ It is dramatically different than the more 
well-studied dynamics of a skyrmion mmm which is 
spontaneously pinned in the absence of external torques. 


Then, Eq. (14) simplifies to 


(av 2 8 \ + V 2 d 2 ) m = mxf + mx (av 2 d 2 — v 2 8 \) m. 

(A4) 


Under the Bogomolnyi relations (Al) the terms on the 


left-hand-side cancel with the last terms on the right- 
hand-side which originate in the damping and the non- 
adiabatic spin torque. Since the same Bogomolnyi rela¬ 
tions are sufficient for the vanishing of the first term on 
the riglit-lrand-side, we conclude that relations (Al) are 
sufficient conditions for solutions of Eq. (A4). 


Therefore, any solution m = n(xi,x 2 ) of Eqs. (Al), 
such as the skyrmion in Eq. (A2), is a traveling solu¬ 
tion m(xi,x 2 ,t) = n(£ i,£ 2 ) of the full Eq. (4j) with spin 
current and damping, with a velocity given by Eq. (A31. 


Regarding calculations presented in Sec. [TV] it is useful 
to note that skyrmion solutions with Q 0 which satisfy 
(Al) have exchange energy IU ex = 4nQ (TUI. Thus, for 


axially symmetric skyrmions, such as (A2), the tensor 

We 


defined in Eq. (121, is diagonal with d^ 


= 5„ 


then see that the velocity in Eq. (A3) coincides with that 


given in Eq. (32) applied for d = 1. 






























7 


[1] A. N. Bogdanov and D. A. Yablonskii, Zh. Eksp. Teor. 
Fiz. 95, 178 (1989). 

[2] A. N. Bogdanov and A. Hubert, JMMM 138, 255 (1994). 

[3] N. Romming, C. Hanneken, M. Menzel, J. E. Bickel, B. 
Wolter, K. von Bergmann, A. Kubetzka, and R. Wiesen- 
danger, Science 341, 636 (2013). 

[4] N. Romming, A. Kubetzka, C. Hanneken, K. von 
Bergmann, and R. Wiesendanger, Phys. Rev. Lett. 114, 
177203 (2015). 

[5] S. Miihlbauer, F. Binz, B. Jonietz, C. Pfleiderer, A. 
Rosch, A. Neubauer, R. Georgii, and P. Boni, Science 
323, 915 (2009). 

[6] X. Z. Yu, Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, 
Y. Matsui, N. Nagaosa, and Y. Tokura, Nature 465, 901 
( 2010 ). 

[7] X. Z. Yu, N. Kanazawa, Y. Onose, K. Kimoto, W. Z. 
Zhang, S. Isiwata, Y. Matsui, and Y. Tokura, Nature 
Materials 10, 106 (2011). 

[8] A. N. Bogdanov and A. Hubert, JMMM 195, 182 (1999). 

[9] A. Leonov, U. K. Rofiler, and M. Mostovoy, EPJ Web of 
Conferences 75, 05002 (2013). 

[10] D. A. Allwood, G. Xiong, C. C. Faulkner, D. Atkinson, 
D. Petit, , and R. P. Cowburn, Science 309, 1688 (2005). 

[11] K. Everschor, M. Garst, B. Binz, F. Jonietz, S. 
Miihlbauer, C. Pfleiderer, and A. Rosch, Phys. Rev. B 
86 , 054432 (2012). 

[12] A. Fert, V. Cros, and J. Sampaio, Nature Nanotechnol¬ 
ogy 8, 152 (2013). 

[13] J. Sampaio, C. Vincent, S. Rohart, A. Thiaville, and A. 
Fert, Nature Nanotechnology 8, 839 (2013). 

[14] N. Nagaosa and Y. Tokura, Nature Nanotechnology 8, 
899 (2013). 

[15] J. Iwasaki, M. Mochizuki, and N. Nagaosa, Nature Com¬ 


munications 4, 1463 (2013). 

[16] J. Iwasaki, M. Mochizuki, and N. Nagaosa, Nature Nan¬ 
otechnology 8, 742 (2013). 

[17] C. Schiitte, J. Iwasaki, A. Rosch, and N. Nagaosa, Phys. 
Rev. B 90, 174434 (2014). 

[18] R. Tomasello, E. Martinez, R. Zivieri, L. Torres, M. Car- 
pentieri, and G. Finocchio, Sci. Rep. 4, 6784 (2014). 

[19] S. Komineas and N. Papanicolaou, Phys. Rev. B 92, 
064412 (2015). 

[20] M. Beg, D. Chernyshenko, M.-A. Bisotti, W. Wang, M. 
Albert, R. L. Stamps, and H. Fangohr, arXiv: 1312.7665 
(2014). 

[21] N. S. Kiselev, A. N. Bogdanov, R. Schafer, and U. K. 
Rofiler, J. Phys. D 44, 392001 (2011). 

[22] J. C. Slonczewski, J. Magn. Magn. Mat. 159, LI (1996). 

[23] L. Berger, Phys. Rev. B 54, 9353 (1996). 

[24] S. Zhang and Z. Li, Phys. Rev. Lett. 93, 127204 (2004). 

[25] M. Finazzi, M. Savoini, A. R. Khorsand, A. Tsukamoto, 
A. Itoh, L. Duo, A. Kirilyuk, T. Rasing, and M. Ezawa, 
Phys. Rev. Lett. 110, 177205 (2013). 

]26] N. Papanicolaou and T. N. Tomaras, Nucl. Phys. B 360, 
425 (1991). 

[27] S. Komineas and N. Papanicolaou, Physica D 99, 81 
(1996). 

[28] J. He, Z. Li, and S. Zhang, Phys. Rev. B 73, 184408 
(2006). 

[29] K. Everschor, M. Garst, R. A. Duine, and A. Rosch, 
Phys. Rev. B 84, 064401 (2011). 

[30] A. A. Belavin and A. M. Polyakov, .JETP Lett. 22, 245 
(1975). 

[31] R. Rajaraman, Solitons and Instantons (North Holland, 
Amsterdam, 1982). 



